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We consider the generating function (prepotential) for Gromov-Witten invariants of 
rational elliptic surface. We apply the local mirror principle to calculate the prepotential 
and prove a certain recursion relation, holomorphic anomaly equation, for genus and 1. 
We propose the holomorphic anomaly equation for all genera and apply it to determine 
higher genus Gromov-Witten invariants and also the BPS states on the surface. General- 
izing Gottsche's formula for the Hilbert scheme of g points on a surface, we find precise 
agreement of our results with the proposal recently made by Gopakumar and Vafa[ll]. 
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1. Introduction 

Since the pioneering work by Candelas et al. in 1991 [6], the theory of the Gromov- 
Witten invariants has been one of the central topics in mathematical physics related to 
string theory. Due to many contributions on this subject we have now well-developed 
mathematical theory [16] [4] of the invariants as well as the concrete methods to calculate 
them applying the mirror symmetry of Calabi-Yau manifolds. However up to very recently 
our concrete methods have been restricted to the genus zero or genus one Gromov-Witten 
invariants. Although we have mathematical definition of the higher genus Gromov-Witten 
invariants, little was known about how to determine them explicitly for a given Calabi-Yau 
manifold. Regarding this a breakthrough has been made recently in [18] for a special class 
of Calabi-Yau manifolds which have a K3 fibration and have a dual description in the 
heterotic string. Independently Gopakumar and Vafa[10][ll] have derived a general form 
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of the prepotential for the higher genus Gromov-Witten invariants, which includes several 
interesting mathematical predictions on the Gromov-Witten invariants. 

In this paper we will propose a recursion relation holomorphic anomaly equation as a 
basic equation for the higher genus Gromov-Witten invariants of rational elliptic surface, 
and will make explicit predictions for them. We find exquisite agreement of our results 
with those by Gopakumar and Vafa. 

To state main results of this paper let us consider a generic rational elliptic surface 
obtained by blowing up nine base points of two generic cubics in P 2 . Under the assumption 
for the cubics the surface S has an elliptic fibration over P 1 with exactly twelve singular 
fibers of Kodaira I\ type. We consider a situation in which the generic rational elliptic 
surface S appears as a divisor in a Calabi-Yau 3-fold X. Since the normal bundle Mx/s is 
given by the canonical bundle K$ we can extract the genus g Gromov-Witten invariants 
N g ((3) of class (3 G H2(S,Z) taking a suitable limit of the prepotential of the Calabi- 
Yau 3-fold X, which is called local mirror principle. Since even for genus zero invariants 
the determination of N g=0 (f3) is technically tedious, in what follows, we will mainly be 
concerned with the following sum of the invariants 

N 9 ;d,n := £ N M 

(f3,H)=d, (P,F)=n 

where H and F represent the pull back of the hyperplane class of P 2 and the fiber class, 
respectively. Associated to these invariants we define generating functions; 

oo oo 
d=0 n=0 

The latter is the genus g prepotential in topological string theory. For g = and g = 1 
we determine it via the local mirror principle applying to X, and find a recursion relation 
satisfied by Z g . n (g = 0, 1; n = 1, 2, • • •) which we generalize for arbitrary g as follows: 
Conjecture 1.1 (Holomorphic anomaly equation for all g) The generating function Z g . n (q) 
has the form 

P 29 +2n-2(0, E 2 , E 4 , E 6 ) (Z 0;1 (q)) n , (1.1) 

with some 'quasi-modular form' for the modular subgroup T(3) of weight 2g + 2n — 2. 
(In the special cases of g = and n = 1, it simplifies to P2 n -2{E 2 , E4, Eq) and 
P2 g (E2(q 3 ), ^(g 3 ), E Q (q 3 )), i.e., exactly the quasi-modular forms of weight 2n — 2 and 
2g, respectively) . And it satisfies the recursion relation 

dZ g . n 1 v^ 1 / w 7 , n ( n + 1 ) v /-, X 

~dE~ = 72 l^ si - n ~ s ) Z 9';sZ g »;n-s + ^ Zg-!. n . (1.2) 

2 g'+g"=g «=i 



We may 'integrate' our holomorphic anomaly equation under certain vanishing conditions. 
In this paper we focus mainly on the special case of n = 1 in which the equation simplifies 
to 



dZ gi i _ 1 
dEo 36 



Zg-i-A ■ (1-3) 
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We integrate (O) with the vanishing conditions and the initial data 



7 f„\ _„§ ggj(gM7) 

^o-Aq)-Q 2 ^3)12 

which has been found in [15] [14]. The Eg theta function comes from the Mordell-Weil 
group of the rational elliptic surface and the eta functions in the denominator come from 
the twelve singular fibers. See [14] for the details and notations. Then we find that the 
solutions Zg-i may be arranged into an all genus partition function of the topological string 
theory: 

Proposition 1.2 (Topological string partition function on S) 



q 



e g8 (3M 7 ) n (l-g 3 -) 4 _ v lfoU 2 ff -2 f2sin A, 2 ( , 



where X represents the string coupling and ti, = e 

We derive the same result following the proposal made in [11] for the BPS state 
counting of the families of genus g curves. From this viewpoint our result ( |1.4j ) comes from 
the following generalization of Gottsche's formula[9] for the Hilbert scheme of g points 
on a surface S: 

Proposition 1.3 (Gottsche's formula with SL(2, C)tX SL(2, C) R Lefschetz action) For 
the Hilbert scheme of g points on a surface S with a fibration structure we can decom- 
pose the Lefschetz SL(2, C) action on H*(S^) into the product SL(2, C)lx SL(2, C) r , 
one for the natural fiber space of and the other for the base space. If we write the 
Poincare polynomial by 

Pt L ,tn(S [9] ) = (tLt R rTr HHSl9 y } (t 2 ^H 2 ^ R ) , 

then the generating function G(t,L,t Rl q) = ^2 g > Pt L ,t R {S^)q 9 , for the surface with 
bi(S) = 0, is given by 



G(t L ,t R ,q) = J] j - 

n>l ^ ^ 



- (t L t R ) n -^qn) (1 - ( tLtR )n+l q n} 

1 , (1-5) 

x 



(1 - V) (1 - ^(^i?) 74 " V) (1 - (t L t R ) n q n ) 



b 2 (S)-2 



We explain our result ( |1.4j ) in terms of the above generalization of Gottsche's formula 

e B8 (3M 7 )G(-*L,-l,f ) =$> fl; i(?)A 2 s- 2 (2sin^) 2 . (1.6) 

This implies that the genus g curves C g in S satisfying (C g ,F) = 1 split into irreducible 
parts, one coming from the Mordell-Weil group and the others from elliptic curves (with 
possible nodal singularities) in the fiber direction. 
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The readers who are not interested in the derivation and the proofs of the holomorphic 
anomaly equation may omit the following two sections and may start from the section 4 
for our main results. 

The organization of this paper is as follows: In section 2, we will introduce a Calabi- 
Yau hypersurface and its mirror, and introduce the hypergeometric series representing the 
prepotential Fo(t) for the Calabi-Yau hypersurface. In section 3, we will take a limit to 
reduce the prepotential to the one relevant to the rational elliptic surface (, local mirror 
principle). We will analyze the reduced prepotential and the mirror maps in detail, and will 
prove the recursion relation, holomorphic anomaly equation at g = 0. Using the formula 
in [1] for Fi, we will also prove the recursion relation at g = 1. In section 4, we will 
propose our recursion relation for all genera, and solve the recursion relation with some 
vanishing conditions. There we also discuss about the Gromov-Witten invariants coming 
from Z g -i(q). In the final section, we discuss some relations to the recent developments on 
the counting problem of the BPS states in topological string theory [11] [15]. There we will 
find a generalization of Gottsche's formula for the Poincare polynomials of S^ 9 '. 
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for her kind hospitality during their stay (August, 1998) at Queen's University where very 
early stage of this work has been done. The research of S.H. and M.-H.S. is supported 
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M.-H.S.), the Ministry of Education, Science and Culture, Japan. The research of A.T. is 
supported by Research Fellowships of Japan Society for the Promotion for Young Scientists. 

After our submission of this paper to the e-print archive, hep-th, we are informed by 
A. Klemm that he is testing the higher genus prepotentials F g for several surfaces other 
than our rational elliptic surface (work to appear). We would like to thank him for sending 
us his preliminary draft prior to publication. 



Note added in proof: In references [15] and [19], a different base F + eg, in stead of our 
H, is used to define Z g= o- n . In this case we consider our generating function Z g - n for the 
invariants N g . d ^ n = N g (/3) summed over (3 with (/?, F + eg) = <i, (/?, F) = n. Then our 
Conjecture 1.1 should be read as follows; 

The generating function Z g . n (q) has the form 

q 2 

P 2 3 +6n-2(£ 2 , E 4 , E 6 ) ^ q y 2n 

with a quasi-modular form P2g+6n-2(-E , 2, E4, Eq) of weight 2g + 6n — 2, and satisfies the 
same holomorphic anomaly equation as 1.2 replaced by the pref actors ^ and 

It is worth while remarking here that in this case the integration constants may be 
determined consistently for all g and n by simply requiring the vanishing conditions for 
the first few terms in the g-expansion of Z g - n (q) : which is defined by 4.5 (, 3.18 and 3.27) 
with D(g, h, k) = C h (g - h, l)k 29 ~ 3 (0 < h < g)[ll}. 
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2. Mirror symmetry of a Calabi-Yau hypersurface in P 2 x Fi 

In this section we will consider a Calabi-Yau hypersurface which contains a generic 
rational elliptic surface as a divisor. We collect necessary formulas for the (genus zero) 
prepotential. 

Let us start with the Hirzebruch surface Fi which is defined by the quotient (C 4 \ 
Z)/ ~ with 

Z = {(ui, u 2 , us, 114) G C 4 I U\ = u 2 = or 113 = 114 = 0} 

and C*-actions 

(tti, u 2 , u s , it 4 ) ~ (Aitti, Ai-u 2 , 1/A1W3, u A ) ~ (ux, u 2 , A2W3, A2W4) , (Ai, A 2 e C*) . 



p2 ~ x 3,75 

X = • 



We may consider a generic hypersurface in the product of the surface Fi with P 2 given 
by the data 

3 

Fi (1,2) 

where (1,2) refers to the homogeneous degrees with respect to the first scaling by Ai and 
the second one by A 2 . The defining equation may be written explicitly as 

^3,3(^1, ^2, ^3,wi,W2)'U 2 + fs,i(z 1 ,z 2 ,zs,u 1 ,u 2 )ul = , (2.1) 

where zi, z 2 , zs represents the homogeneous coordinate of P 2 and #3,3 (, ) refers to a 
generic homogeneous polynomial with bi-degree (3,3) (, and (3,1), respectively,) for the 
coordinates zi,z 2 ,zs and u\,u 2 . This is an elliptic Calabi-Yau hypersurface over Fi with 
the Hodge numbers h ' = 3 and h 2,1 = 75. Two of the three elements in H ' (X) come 
from the base Fi and the other comes from the fiber elliptic curve. We may find in X 
a rational elliptic surface £ with its defining equation of bi-degree (3, 1) in P 2 x P 1 . It 
appears as a divisor u% = 0, which is the cubics in P 2 over the (—1) curve in Fi. 

The positive classes in H 2 (X, Z) are generated by the three integral elements in 
i/ 1 ' 1 (Y) corresponding to the divisors 

H = ( Zl = 0) n x , F = ( Ul = 0) n x , D = (« 4 = o)ni . 

We sometimes denote the corresponding forms by J\ , J 2 and J3 , respectively. It is straight- 
forward to determine the non-zero intersection numbers K^ c = j x J a A A J c and 
c 2 J a = J x c 2 (X) A J a with the second Chern class c 2 (X) to be 

Tftop _ r, Tftop _ o Tftop _ o 17-top _ o 

-^112 ~ ^ 1 -"-113 — ' xv 123 — ' -""133 — ' 

c 2 Ji = 36 , c 2 J 2 = 24 , C2J3 = 36 . 

The ambient space P 2 x Fi is so-called the toric Fano manifold, and thus we can easily 
construct the mirror Calabi-Yau hypersurface X y based on Batyrev's toric method[Bat]. 
Furthermore the prepotential of the mirror Calabi-Yau hypersurface X y is determined by 
the general formula obtained in [12] [13]. 
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Here we collect necessary formulas to determine the prepotential. We start with 
a hypergeometric series representing a period integral for a deformation family of X v 
parameterized locally by x, y, z; 



wq(x) = c(n,m,k)x n y rn z k 



with 

T(l + 3n + m + 2k) 

c(n, m, k) — 



r(i + n) 3 r(l + m) 2 r(l + k - m)r(l + k) 

The local parameters (x,y,z) has been chosen so that its origin represents the celebrated 
boundary point where the monodromy is maximally degenerated [21] [22]. The series wq(x) 
represents the period integral for the invariant cycle about this degeneration point and sat- 
isfies Picard-Fuchs differential equation (see Appendix). As a complete set of the solutions 

of the Picard-Fuchs equation, we have {wq(x) , Wa (x) , (x) , w^(x)} (a, b = 1,2,3) 
where 



with u>o(x, p) = J2 n m k>o c ( n + Pi,™ + Pit k + p 3 )x n+Pl y m+P2 z k+P3 . In terms of the 
solutions of the Picard-Fuchs equation, the mirror map is defined by the relation 

t a = ^j^- (a = 1,2, 3), (2.2) 

which connects the deformation parameters (x, y, z) to those (ti, £2, H) parameterizing the 
complexified Kahler moduli of X at the large radius (Im(£ a ) — > oo). Now the prepotential 
of the mirror X v is defined to be 

F{2) = { w " (wm ~ ? ^ w > )] + ? 

Then the mirror symmetry conjecture asserts that the prepotential F(x) of the mirror X v 
combined with the mirror map provides, up to terms of classical topological invariants, the 
generating function Fo(t) of the Gromov-Witten invariants of X; 

Fit) = - V K top t Ut -V ( c 2 J fa)^ , C(3)x(A) 1 N (8)e 2ir ^T, W 

[ ) 3\ 2^ ^abcWc 2_. 2 4 + 2 (2ttz)3 + (2ttz)3 2^ JV ^ e 

a,b,c v ' v ' 0^/3eH 2 (X,Z) 

where we substitute the inverse relation of ( j2.2| ) into F{x). 
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3. Holomorphic anomaly equations for g = 0, 1 



3.1. Local mirror principle and the reduction to S 
Let us consider the following limit 

F(q,p) := F(t) - (topological terras)|i mt3 _K» 

with q = e 27 ™* 1 and p = e 27 ™* 2 . Since the class J3 measures the volume of the fiber 
P 1 of Fi parameterized by ^3,^.4 and the volume of the curve contained in the divisor 
S = (W3 = 0)nJ are measured to be zero by this class, the limit Inr^ — > 00 throw away all 
the Gromov-Witten invariants except those of the curves contained in the rational elliptic 
surface S. Thus we may expect that the reduced prepotential F(q,p) coincides with the 
generating function defined in Section 1. This is so-called the local mirror principle, and 
somehow generalize the arguments done for the isolated (—1, — l)-curves in Calabi-Yau 
manifolds [17] [8]. 

In our case of the elliptic Calabi-Yau hypersurface X in P 2 x Fj, the limit Imt3 — ► 00 
which translates to the limit z — > in the mirror X y greatly simplifies the period integral 

w (x); 

. s-^ T(l + 3n) n 

w (x, y, z)\ z=0 = 2^ rri )3 X ' 

n>0 1 ' 

=: <j>(x) . 

We note that the series 4>(x) is nothing but a solution of the Picard-Fuchs equation of the 
fiber elliptic curve of S: 

{9l - 3x(3# x + 2){W X + l)}0(x) = . (3.1) 
Another solution of (|3.1| ) about x = may be given by 

fa) = \og{x)${x) + pS! Z ^3 (^(1 + 3n) - 3^(1 +n))x n , 

^ r(i + n)3 

where ip(z) = -4-logT(z). 

Now looking at the relations between the 'periods' of X via the prepotential and the 
period integrals of X v \ 



it is straightforword, although involved technically, to derive the following consice form for 

d 
dU 



the derivative -^-F(p, q) (t v := i?) under the limit Imt 3 — > 00 (cf. Sect. 7 of [14]). 



Proposition 3.1. 

— F(q p) = T^V n 
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where we define 

f n (x) = --{4>(x)£ n (p(x) - <fi(x)£ n (p(x)} , 

with a linear operator 

(_-\\n n 

77, x n! - LJ - 

k=i 

Remark 3.2. By constructing Barnes integral representation of the series 4>(x), it is an 
easy exercise to make an integral symplectic basis of the Picard-Fuchs equation (|3.1|) . It 
turns out that our bases (j)(x) and <p(x) in fact constitute an integral symplectic basis 
about x = 0. Therefore we may write the holomorphic one form of the elliptic curve by 
= <PqA + (piB, where A and B are symplectic bases of the elliptic curve. Then the 
function fi(x) may be written by —3 J e ^e A 9 x Qe- This is so-called the classical Yukawa 
coupling [6] of the elliptic curve, and may be determined from the Picard-Fuchs equation 
( j3Tl) to be fi(x) = 1 _ 9 2 7 X - -^ or ^ ne other functions fm{x){m > 2) we will find a powerful 
recursion relation. 

The relation ( |2.2| ) is also simplified in the limit z —>■ due to the following relations 

w o (x,y,0) = (f)(x) , w^ix^y.O) = (j)(x) , w£\x, y, 0) = + C m (f)(x)y rn , 

m>l 

where £(x) = ^2 n>0 j^rjl (ip(3n + 1) — "0(1)) x n . As the inverse relations of ( |2.2|) , we have 
a; = x(q,p) and ?/ = y(q,p) with g = e 771 u, o(*.w.o) and p = e 7ri u, o(*.w.o) . 



Proposition 3.3. Under the limit z — > 0, we /md £/ia£: 

(%) TTie inverse series x(q,p) does not depend on p and is given by the level three modular 
function; 

x{q ' p) = • (3 - 2) 

12/ \ 

where Hb{q) = Ji2( g 3) ^ s ^ e Thompson series in the notation of [7]. 
(2) The inverse series y(q : p) is determined iteratively as a power series of p through the 
relation 

y(q,p) = pip(x)e S m >i Cm ( x ) y ^ (3.3) 



where ip{ x ) = e * (x) an d c m(x) = C ^}^ 



Remark 3.4. The function ip(x) with x = t3 ^ q \ has first appeared in [14], and has been 
determined in terms of the modular functions of level three[28]; 



2 
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where £3^4(9) = x (i-27x) • Also the following relations are standard results coming from 
the Gauss-Schwarz theory for the Picard-Fuchs equation flOl); 

<P(x(q)) = e 3 (q)8 3 (q 3 ) + 8 2 (q)8 2 (q 3 ) , (3.4) 

i (l 1* 

<P{x) 12 x{l-27x) 3 = 77(g) 24 , = </>(x) 2 x(l - 27x) , (3.5) 

where 8 2 (q) = E m6Z 9 (m+ ^ )2 and 3 (g) = Emez^- 



The following lemma may be derived directly from the definition c m (x) and the relation 

8 x <j>(x) 1 / _ h fi E 2 (q) 



4>(x) 3 V 12 36</>(x) 2 

which follows from ( |3.5| ) ; 

Lemma 3.5. Under the relation x = t3 ^ q - ) , the function c m {x) = Cr ^^ fnay be written 
by 

c m (x) = B m {h)^% + D m {h) , (3.6) 

where B m and D m are some polynomials of f± determined by the following recursion rela- 
tion; 



D m+ x = ~ (m 7i)2 - 8 ) S - + ^ + m + ^} D ^} 



(3.7) 



t/jit/i initial values B\ = — 1| and Di = — 1|. 



We present here the first few terms coming from the recursion relation ( |3.7| ); 

B2 = ^ A 2 , b 3 = ^ (a 2 - ^ a 3 ) , b, = (a 2 - I a 3 + J A 4 ) , • • • (3.i 

1)2 = "3^ ( Jl ~ 4O ' ^ 3 = ~ ik ( A ~ + ^ 

D = JL />_?/? 3/T (J<)) 

4 3888 V 17 68 
We can verify directly using the Picard-Fuchs equation ( |3.1| ) that the formal solution of 
the recursion relation may be written in terms of the functions f m (x) in Proposition 3.1 

as 

B m = ~^ , D m = i- ( (m + 1)2 f m+ i + (3g, + m + 2 - ■ (3.10) 

36 ji [ m 12 J 



As a result we see that the functions / m (x)'s may be determined in terms of the recursion 
relation ( p.7|) . 

Since both the B m and D m are polynomials of fiix) = 1 _ 9 2 7 X = , they have nice 

behavior under the level three modular subgroup T(3). Therefore the modular anomaly 
comes from the E^-term in c m . We may express this anomalous behavior via the partial 
derivative of c m ; 

dc m (x) ^ 1 fm{x) , v 

8E 2 36</>(x) 2 ' 1 ' 

which plays a central role in the following derivations of the holomorphic anomaly equa- 
tions. 

3.2. Holomorphic anomaly equation at g = 

Now we are ready to prove the recursion relation for Zo ;n (q , )'s, which come from the 
mirror symmetry conjecture through the expansion 

^.rt = E^" = E»W». (3.12, 

p n>l rV ' n>l 



Theorem 3.6. (Holomorphic anomaly equation at g = (c.f. [19])) The function Zo- n 
satisfies the recursion relation 

dZ 1 v— -i 

= ^5^( n - S ) sZ 0;n- S ^0; S (fl > 1) . (3.13) 

2 s=l 

(Proof) From Proposition 3.1 and Proposition 3.3, we have 

where the quasi-modular property (anomalous modular property) appears in c m through 
( [3. 11|) . Now we first note 

d y 

-y 



dt p i + Em=i mc ™y 

Using this and the relation (|3.11|) , we have 

dy -yY2=i%my m 1,0 w s 



^2 l + Em=i TOC »# m ZQ x dt/> K dt 
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(is-vXoHO • ( 3 - 14 ) 



Now we have 

L - x f m {x)my 



m>l 



m>l 



which says, up to constant terms for p, that 

F(q,p) 



dF(q,p) 1 f d 



<9£ 2 72 
This proves the recursion relation. 



Now we determine the explicit forms Zo- n for lower n from the formula ( |3.12| ). After 
coming to a conjecture about the form of Zo- n , we remark that the holomorphic anomaly 
equation ( |3.13| ) with some obvious inputs for the Gromov-Witten invariants suffices to 
determine the form Zo- n for all n. 

Since our formula ( |3.12| ) is written in terms of the known functions fi(x), ip(%) with 
x = tsB / q \ and E2(q) for each order of p, and also the order of p coincides with the order 
of ip(x), it is easy to deduce that Zo- n has the following form in general, 

Z Q ;n = Co ; n(/l,§)0 2n ( , (3.15) 



where we have factored the form of Zq-i = Looking into the detail of the expansion 
of ( ^.12| ), we see in general that <fi 2 x Go- n (fi, -p-) is a polynomial of -|f with coefficients 
being polynomials of l//i over Q. Here we present the first few of them, 



Gq-2 = 






- 1 f 


31104 \ 


Gq-5 = 


269 


62208000 



Go . 3 = JL _ A)^ + 3 2 

7776 V h 5 2 

2000 12 24 2 3125 A _ 10 
~2m {1 ~J + J! ] E2<t> + 21789^ 
Now we note the following relations for the polynomials of l//i; 

E A =9 (8- -^V , E G = -27 (l - ^ + ^ . 

V hJ \ h ft) (3.16) 

270 8 - 18£ 4 </> 4 -E\- 8E 6 <p 2 = 
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Using these relations we find that 



Gq-s =—^- \E 4 + —EV\ d)- 6 , G -4 = — - — (E 6 + 5E 4 E 2 + 12E%) </r 8 , 
' 69984 V 5 2 J ' 839808 V 2> r 

G .5 = ( —E\ + -E 6 E 2 + —E4E0 + 25El) 0~ 10 . 

0,5 40310784 V 125 4 3 3 2 2 J 

Here we observe explicitly that disappears nontrivially in the final expression of Zo- n for 
lower n's. We do not have general proof about this but may state it as follows; 



Conjecture 3.7. The function Zo- n (q) in (W^TE) takes the form 



Z . )n (q)=P 2 n-2(E2,E 4 ,Ee)(Zo.i{q)) n , (3.17) 
where P 2 n-2 is a quasi-modular form of weight 2n — 2. 



Remark 3.8. The function Zo- n contains the multiple cover contributions. We may 
subtract these contributions considering 

Z ., n (q) := Z . n (q) - £ Lz , n/k (q k ) . (3.18) 

k\n,k^l 

The q series coefficients of Zo- n "count" the numbers of rational curves C in our rational 
elliptic surface S satisfying (C, H) = d and (C, F) = n. The homology classes of curves in 
S, in general, have the form [C] — dH — a\e\ — a 2 e 2 — ■ ■ ■ — ageg with ai, a 2 , ■ ■ • , ag > 0. 
Therefore we have 3d — a± — a 2 — ■ ■ ■ — a$ = n for (C, F) = n, which implies 3d > n. In 
other words, we should have 

~ n 

N ., d;n = for d< - (3.19) 

for the coefficients of Zo- n . From a simple counting of the dimensionality of the quasi- 
modular forms of weight 2n — 2, we see that the vanishing condition Q3.19|) together with 
the above Conjecture 3.7 provides sufficient data to determine the integration constants 
for the recursion ( |3.13j ), and determine completely Zo ;n for all n (, see Section 4). 



3.3. Holomorphic anomaly equation at g = 1 

According to [1], we have an explicit expression for the genus one prepotential 
Ff cov (t) in terms of the discriminant of the hypersurface X; 

F* cov {t) = logj (^^J cfe(*,y,z)-^-V 3 z- 4 det } , (3.20) 
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where the discriminant may be determined from the characteristic variety of the Picard- 



Fuchs equation presented in Appendix. Several exponents in (|3.20|) have been fixed by 
the requirements of the asymptotics of F BCOV when Imt a — > oo. The explicit form of the 
discriminant is a complicated polynomial of x,y and z, however, in the limit z — > 0, it 
simplifies to 

dis(x,y, z)\ z=0 = (1 - 27x) 3 {l - 27x + (1 + yf - 1} . (3.21) 

Also under this limit, it is easy to show from the definition that the mirror map 

z(q,p, r) (r = e 27 ™* 3 ) simplifies to 



2g(aQ | c m (x)y r 

z(q,p,r)=re ^-"-»>i m 

= rip(x) 



^3 P 



(3.22) 



y(q,p) 



Now using the relations (|3.5|) it is straightforward to derive; 



Proposition 3.9. When r — > 0, up to the topological term — J2 c ( c ^^c)t c = — 31ogg — 
2 logp — 31ogr, we have; 

Fr° V {q,p) =log {0 4 (1 - 27s) M V (q)-^} 

L e ~E m>1 ^y m (3.23) 
+ log {{(1 - 27.) + (1 + yf - _ } . 



Remark 3.10. 

1) For the p-independent term of F BCOV (q,p) the local mirror symmetry does not apply. 
This is because these curves are parallel to the fiber elliptic curves of X and therefore 
can move outside of the rational elliptic surface. In fact we see that 



N, 



s=l;d,0 



4 d = 3 
d ^ 3 



after subtracting the genus zero contributions N g =o-d,o = 168 (d = 1,2 mod 3), 144 (d = 
mod 3). The number 4 should be regarded as the Euler number of the base Fi for 
the elliptic fibration. 

2) There is a difference in the normalization of the prepotentials between F BCOV and 
our F g =i in the introduction. These are related by the factor 2 coming from the 
orientation of curves as 

F 1 (q 1 p) = ±F 1 BCOV (q 1 p) . (3.24) 



Now we define the generating function Zi- n through 

F 1 (q,p) = J2 Z iAQ)P r 
p>i 
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and prove the holomorphic anomaly equation at g = 1. 



Theorem 3.11. (Holomorphic anomaly equation at g = 1) The function Zi- n (q) satisfies 
the following recursion relation; 

dZ^n 1 ^ n(n + l) 

~QE~ = 36 l^( n ~ s > s Zl ' n ~ s Z °' s + 72 0;n ^ n - 1 )' ^ 3 - 25 ) 

2 s=l 

(Proof) Since we have already shown the relations 

dy 1,9 w 9 9c m 1 / m 



9£ 2 36 v 9t/ /v 9t p u/ ' 9£ 2 36</> 2 ' 
and 

9 = y 

9t p 1 + Em=l mC mV m 

it is straightforward to derive 

9^2 9y 1 + L m >i TOC m(£)y 9-E 2 



: <W (l + £ m >i™c m2 /™) 

^ m>l 

1 / 9 e -E m >i c -( x ^ m \ a 

= 4 Nr log ft 1 " 27x ) + ( X + ^ " ^"' TTv TT^} x 

36 I 9t p 1 + 22 m >i mc ™( x )y I dt P 

+ 36 9^ Fo + 36 ( 9^ ) F ° 
19 9 19 9 

= 36 ( 9^ FiBCW)( 9i; Fo) + 36 9i; ( 9i; + 1)Fo 

Taking into account the difference of the normalization ( |3.24j ), we conclude the recursion 
relation. I 



Now we may determine the generating function Z\- n (q) explicitly from ( |3.23|) under 
the relation Fi(q,p) = hFf cov (q,p). As in the cese of genus zero, we may represent Zi ]n 
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in terms of fi(x), ip{x) and E 2 (q). Corresponding to (|3.15|) , we have 



Z 1;n = G ls „(/i, ||)</> 2n (f^Y = G 1 , n (fu |f )</> 2n (Zo-^r . (3.26) 



After straightforward evaluation of (|3.23 ), we obtain for the first few of Gi ;n 's; 



' 15552 VV A 2 /^ 6 V 13/i/ ^ 2V ^ 18 V 
If we use the relations ( |3.16| ) for the polynomials of l//i , we find 

Gl;l = ^ ( # 2 + ^ #" 2 , G li2 = 77^ (V - i«4 + ^ 2 -E 2 + 1 ' 



S + 5*) • G - =452 (** " + I ^ + B- 



3 \ . 



Contrary to the case of (7 = 0, the ^dependence remains in Zi- n (q) after the elimination of 
the polynomials of l//i. Thus we arrive at the follwing weaker statements about Zi- n (q); 



Proposition 3.12. The generation function Zi jTl (q) in ( \3. 26]) takes the form 



Z 1 , n (q) = P 2n ( ( j )l E 2l E 4l E 6 )(Z . 1 (q)) n , 
where P 2n is a 'quasi-modular form' of weight 2n for the modular subgroup T(3). 

Remark 3.13. 

1) The form of the polynomial P 2n is not unique because of the relation ( |3.16| ) among 
<fi, E4 and Eq. 

2) The function Zi- n contains the contribution from the genus zero curves (i.e., the 
degenerated instanton [1][8],) as well as the contribution from the multiple covers. We 
may separate these in Zi- n as follows; 

Z 1 . n (q) = Z 1 . n (q) + ^-l(k)Z 1 . n/k (q k ) + ^J2^Z ., n/k (q k ) , (3.27) 

k\n,k^l k\n 

where cr_ 1 (A;) = J2 m \k m - ^ ne function Zi- n (q) is expected to 'count' the numbers of 
the elliptic curves C with (C, F) = n in S. As is the case of g = 0, we have certain van- 
ishing conditions for the elliptic curves, which is useful to determine the integration 
constants for our holomorphic anomaly equation ( 3.25|) . However as we will argue in 



the next section the appearance (j) in the polynomial P 2n increases unknown param- 
eters in the integration constants. From this reason holomorphic anomaly equation 
become less powerful than the case of g = to determine Z\. n . 
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4. Predictions for Gromow-Witten invariants of higher genera 

4-1. General considerations 

From the analysis for g = and g = 1, we naturally come to the following conjecture 
about the holomorphic anomaly equation for all genera. 



Conjecture 4.1. (Holomorphic anomaly equation for all g) The generating function 
Z g . n (q) has the form 

P 2g+ 2n-2(^E 2l E 4l E 6 )(Z . 1 (q)) n (4.1) 

with some 'quasi-modular form' for T(3) of weight 2g + 2n — 2. (In the special case of 
g = 0, it simplifies to P 2 n-2(-E , 2, E±, E G ), i.e., exactly the quasi-modular form of weight 
2n — 2). And it satisfies the recursion relation 

dZ g . n 1 ^ ^ n(n + l) 

~dE 2 ~ = T2 ^ }_^s{n- s)Z gl , s Z gll . n _ s + Z g _ 1;n . (4.2) 

g'+g"=g s=1 



In the following we will consider the solutions of the holomorphic anomaly equation. 
For this purpose, first of all, we need to have some data to fix the 'integration constants' 
for our recursion relation ( |4.2| ). Let us suppose that a curve C g in the rational elliptic 
surface is in a homology class [C g ] = dH — aiei — a 2 e2 ■ ■ ■ — ageg, where e^'s refer to the 
— 1 curves from the blowing ups. Then, since Z g - n counts the Gromov-Witten invariants 
of genus g curves with ([C g ], F) = n, we should have 



and for the arithmetic genus 



([C g },F) = 3d-J2 



(1; 



n 



(4.3) 



9a(Cg 



(d-l)(d-2) 



9 

E 

1=1 



a>i(a>i 



(4.4) 



If the curve C g is singular, the arithmetic genus might be different from the genus of 
the normalization of C g . We will come to this point later, however for the moment we 
will ignore this difference. Then it is easy to see that the above two constrains pro- 
vide us several vanishing conditions on the numbers of curves. In Table 1 we have pre- 
sented the lowest degree d = ([C g ],H) for which a curve C g may exist for given g and n. 



n\g 





1 


2 


3 


4 


5 


6 


7 


8 


9 


1 





3 


6 


9 


12 


15 


18 


21 


24 


27 


2 


1 


3 


4 


6 


7 


9 


10 


12 


13 


15 


3 


1 


3 


4 


4 


6 


7 


7 


9 


10 


10 


4 


2 


3 


4 


4 


5 


6 


7 


7 


8 


9 


5 


2 


3 


4 


4 


5 


5 


6 


7 


7 


8 



Table 1. Each number shows the lowest degree d 
given g and n. 



([Cg],H) for curves of 
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The homology classes of curves in the table may be written explicitly for given g and 
n. For example, for n = 1 they are simply given by [C g ] = + gF (i = 1, • • • , 9) . The 
data in the table provides us vanishing conditions for the first few terms in the g-expansion 
of Z g - n (q), where tilde represents the subtraction of the degenerated instantons from the 
Gromov-Witten invariants. We may relate Z g - n (q) (g > 2) to the subtracted functions by 

Z g . n (q) = Z g , n {q) + D(g,g,k)Z g . n/k (q k ) + J2z2 D ^> h i k ^h-,n/k(q k ) , (4-5) 

k\n,k^l h=0 k\n 

with some rational numbers D{g 1 h, k). Therefore all we need to fix from the vanishing 
conditions are the 'integration constants' together with the rational numbers D(g,h,k). 
For higher n > 4 it turns out that the vanishing conditions in Table 1 are not sufficient to 
determine completely both the integration constants of the recursion (|4.2[) and the form 
of the degenerated instanton D(g, h, k). However we will see based on a simple counting 
arguments that for lower n < 3 they suffices at least to fix the integration constants for 
our recursion relation. Especially for the case n = 1 they determine both the integration 
constants and the form of the degenerated instanton. 

As an extreme case let us first consider the genus zero generating functions 
Zo. n (q) (n = 1,2, •••). As we have already considered in Remark 3.8, the only unde- 
termined in this case are the integration constants in the polynomial -P2n-2 in ( |3.17| ). We 
see that these constants can be fixed by the vanishing conditions in Table 1 comparing the 
first column of the table with the dimensionality of the modular form of a given weight, 
which is given by the series 

1+t 4 + t 6 +f 8 +t 10 + 2tl 2 +t 14 + 2t 16 + 2t 18 + 2t 20 + ... (46) 



(l-t4)(l_ £ 6) 

Now let us look at the functions Z g]n (q) for all g and lower n(< 3). In this case we only 
have the weaker assumption on P2 g +2n-2 i n Q4.1|) . Taking into account the relation among 
<p, £4 and Eq in the second line of ( |3.16| ), we may estimate the relevant dimensionality 
of the modular form P2 g +2n-2\E 2 =o G integration constants of the holomorphic anomaly 
equation,) by 

( n + - 2 4 tn 4 - + 6? =l+t 2 +2t 4 + 3t 6 + 3t 8 +4t 10 + 5t 12 + 5t 14 + 6t 16 + 7t 18 + 7t 20 + --- . 
(1 t ){1 t ) 

(4-7) 

The growth of the dimensions should be compared with each line of the table 1 under a 
suitable shift. From these comparisons of the numbers of the 'integration constants' and 
the numbers of the vanishing conditions, we may deduce that for n < 3 the vanishing 
conditions suffices to determine the integration constants while leaving some of D(g, h, k) 
undetermined. To go beyond this rather unsatisfactory situation, we need to know more 
details about the numbers of genus g curves of a given homology classes or the form of the 
degenerated instanton D(g, h, k). Though our simple vanishing conditions are insufficient 
to determine all the unknowns, we see from the first line of the table they are restrictive 
enough to fix the form of Z g .i(q) and D(g, h, 1) completely. In the next subsection we will 
present a detailed analysis of Z g -i(q) for all g. 
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4-2. Gromov-Witten invariants Z g -i and degenerated instantons 



To get some intuition about the curves C g in S, let us first recall the form of Zo ; i(g) 
obtained in [14] [28] 



where 7 = (1, - - - , 1, — 1) . The appearance of the E% theta function originates from the 
well-established fact that the sections of the rational elliptic surface form additive group 
called Mordell-Weil group, and it becomes a lattice isomorphic to the E% lattice endowed 
with a positive definite bilinear form[25] [14] [24]. The eta functions in the denominator 
have been explained by introducing pseudo-section which is a composite of a section with 
some of the twelve singular fibers with its homology class a + kF. Thus the function Zq-i 
counts the numbers of pseudo-sections in S which are not irreducible but naturally come 
in the theory of the stable maps[5]. 

Now in our general case of genus g curves with (C g ,F) = 1 the function Z g -i counts 
the genus g sections of the elliptic fibration of S. Since the generic fiber spaces are elliptic 
curves, the genus g section are composite of two components, one is a pseudo-section and 
the other consists of g fiber elliptic curves. The genus g sections of the lowest degree are 
those with their homology classes given by 



[C g } = e i + gF (i = l,.-.,9). (4.9) 

Thus the expansion of Z g -\(q) start from q 3g with its coefficient 'counting' the number of 
the genus g sections of class ( [Of ) . (This is the vanishing condition we have listed in Table 
1 for n = 1 and g.) The g elliptic curves in general avoid the twelve singular fibers and 
make a ^-dimensional family parameterized by Sym 9 (P 1 ). 

As we have already remarked, the vanishing conditions for the 'numbers' of curves grow 
much faster than the dimensionality of the integration constants ( [4.7| ) plus the numbers of 
the unknowns D(g, h, 1) in 



Z g . 1 (q) = Z g . 1 (q) + J2D(g,h,l)Z h>1 (q) . (4.10) 

h=0 



Owing to this nice property, we can integrate our holomorphic anomaly equation for n = 1, 

:Z g -i;i , (4.11) 



dZ g , x _ 1 



dE 2 36 

with the results listed in Table 2 for Z g -i(qys. For the degenerated instantons D(g, /i, 1) 
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we find 





= z 1;1 


12 




^2-1 


= Z%;\ 


+ X (M 2 )Z . 1 




^3-1 


= Z 3 -i 


+ ^x(M 3 )Z ;i - 


12 ' 


Z^-l 


= Z4.1 


+ i X (M 4 )Z 0;1 + 


1 Z 

360 ' 


Zb-,1 


= Z 5 -i 


+ ^x(M 5 )Z 0;1 - 


1 z 

20160 



(4.12) 



6 " 



1 * U 

1 — 7^4,1 



80 ' 

with the orbifold Euler number of the moduli space of genus g stable curves, x(M ff ) = 

I B I 

2 g (2g-2)! (d — 2)- In fact these forms of the degenerated instantons -D(g, /i, 1) (0 < /i < gf) 
coincide with recently established results in [11] [23], where we have 

9 

Z g ., 1 (q) = Y, C h(g-h,l)Z h;1 (q) , (4.13) 

h=0 

with the coefficients determined by 



Zq-i — 


9 + 36q + 126q 2 + 360q 3 + 945g 4 + 2268q 5 + 5166g 6 + 11160g 7 + • • • 


z 1;1 = 


-18g 3 - 72q 4 - 252q 5 - 774q 6 - 2106g 7 - 5292g 8 - 12564g 9 


Z2;l = 


27q 6 + 108q 7 + 378q 8 + 1224q 9 + 3411q 10 + 8820g n + 21663q 12 + • • • 


Z3-l = 


-36q 9 - U4q 10 - 504q 11 - 1710g 12 - 4860? 13 - 12852g 14 


Z±-l = 


45q 12 + 180q 13 + 630q 14 + 2232g 15 + 6453g 16 + 17388q 17 + • - • 


Zb;l = 


-54q 15 - 216q 16 - 756q 17 - 2790g 18 - 8190g 19 - 22428g 20 



Table 2. Solutions of the holomorphic anomaly equation (|4.11| ). These are 
related to Z g;1 by (P^ ). 



We note in Table 2 that for the first three terms in the expansion we have 

Z g ,i{q) = (-l) 3 X (Sym 3 (P 1 ))(9^ + 3§q 9+l + 126q 9+2 + ■ ■ •) (4.15) 

where x(Sym s (P 1 )) represents the Euler number of Sym s (P 1 ) = P 9 . This is also in 
agreement with the argument in [11] for counting curves with moduli. From the fourth 
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term in the series expansion ( 4.15 ) contributions from the singular fibers come in, which 
somehow generalize the situation we encountered in the case of g = 0. 

For the first few of Z g -x (q) 's we have determined explicitly the forms of the polynomials 
P2g((f>, E 2 , E4, Eg). We find that if we use the following relations 



3E 2 (q 3 ) =2<P(q) 2 + E 2 (q) , 9E 4 (q 3 ) = lO^q) 4 - E±(q) 
27E 6 (q 3 ) = 35</>(?) 6 - 7£ 4 (^(<7) 2 ~ E 6 (q) , 



(4.16) 



they summarize into concise forms; 



Q 2e E8 (3t,t 7 ) 1 3 

z o-,i - — ^(^3)12 — ' - 12 J 0;1 

^i = I ^ 5 (5^(g 3 ) 2 + ^( ? 3 ))^o ; i ^ 4 - 17 ) 

^3;i = ^^{35E 2 (q 3 ) 3 + 21E 2 (q 3 )E 4 (q 3 ) +4E 6 (q 3 ))Z . 1 . 

Proposition 4.2. The solutions of the holomorphic anomaly equation (\4-U\ ) ^ a ^ e the 
following general form 

Zg-M = P2 9 (E 2 (q 3 ), E 4 (q 3 ), E 6 (q 3 )) Z Q , x {q) , 

where P 2g is a quasi-modular form of weight 2g. 

The reason of this simplification will be explained in the next section. Here for later 
use we define Q g; i by 

Z g -i(q) =e E8 (3t,t 1 )g g . 1 (q 3 ) , (4.18) 

which should count the genus g pseudo-sections made from a section, say the zero section. 
We verify directly that the functions {7 5; i's depend on q through q 3 and have the following 
expansions; 

£o ; i = 1 + 12 q 3 + 90 q 6 + 520 q 9 + 2535 q 12 + 10908 q 15 + 42614 q 18 + ■ ■ ■ 

g vl = - 2 q 3 - 30 q G - 260 q 9 - 1690 q 12 - 9090 q 15 - 42614 q 18 

g 2 . 1= 3 q 6 + 52 q 9 + 507 q 12 + 3636 q 15 + 21307 q 18 + 107772 q 21 + ■ ■ ■ 

g 3;1 = - 4 q 9 - 78 q 12 - 840 q 15 - 6570 q 18 - 41580 q 21 - 225432 q 24 

g 4;1 = 5 q 12 + 108 q 15 + 1271 q 18 + 10756 q 21 + 73083 q 24 + ■ ■ ■ 
£ 5;1 = - 6 q 15 - 142 q 18 - 1812 q 21 - 16494 q 24 - 119770 q 27 

In the next section we will discuss geometric interpretation of the numbers in the above 
expansions. 



(4.19) 
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5. Discussions 



5.1. Counting BPS states 

In this section we would like to discuss relations of our results to the very interesting 
proposals made in the recent works by Gopakumar and Vafa [10] [11]. 

It is known in physics that the genus g prepotential F g has a meaning as the genus g 
topological partition function of the twisted Calabi-Yau sigma model in the type IIA string 
theory and it does not receive the string perturbative corrections due to the fact that the 
dilaton belongs to the hypermultiplet in the type IIA string. Since the heterotic/type II 
string duality connects the heterotic dilaton field to one of the vectormultiplet moduli in 
the type IIA side, we may expect to extract the non-perturbative properties in the type 
IIA side from the perturbation theory in the heterotic string. In our case of the topological 
amplitude F g , Gopakumar and Vafa [10] [11] have found that it can be derived from the one 
loop integral in the heterotic side. They found that the Schwinger one-loop calculation for 
a particle in a constant background electro-magnetic filed applied to the BPS states with 
spin (ji, J2) under SO {A) = SU(2)l x SU(2)r determines the higher genus F g . According 
to [11] the contribution of each BPS state with spin (j'1,,7'2) to F g is determined by the 
following decomposition with respect to SU(2)l C SU(2)l X SU(2)r; 

2ji 

(ji,j 2 )=l> r [(-)©2(0)r , (5-1) 

r=0 

where we allow a r formally to be negative integer. Then the topological partition function 
has the following expression 

00 2 il( r ) r oo 7 2r-2 

E A25_2 ^= E EE^( r )/ -( 2sin |) e-^Cn-Hm) , (5.2) 

g=0 F-.BPS r=0 m€Z " e S 

where A measures the central charge of the BPS state F and the summation over m is 
explained as the central charge of the fifth dimension which originates in the M-theory. 
In the type IIA theory the central charge A is measured by the Kahler classes of the 
corresponding Calabi-Yau manifolds, and after the integration we arrive at 

00 1 / h \ \ 2r_2 

E^w= E E5>4( 2sin Tj qk " • (5J) 

5=0 0^ v eH 2 (X,Z) r>0 k>0 V ' 

where ofy represents the summation of a r (T) over the BPS states V with charge rj, and 
qkij _ exp(— 27r(/c?7, K)) with K = t\ J± + • — h t h i,\J h \,\ the Kahler class. Expanding ( |5~3| ) 
with respect to A we obtain 

FM = E E Y, a H k k 29 -"C h {g - h, , (5.4) 

0^ v <EH 2 (X,Z) h=0 k\ v 
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which is the general form proposed in [11] and we have reproduced in our special case 
see and (PI)). 



In the type IIA picture the BPS state appears as the D2 brane with the flat U(l) 
connection. Let us suppose that a genus g curve in a Calabi-Yau manifold comes with 
the moduli of the deformation M.. Together with the Jacobian of each curve we have 
the fiber space M. — > M. for the BPS states. Gopakumar and Vafa argued based on the 
M-theory that in this case the counting BPS states in ( |5.2| ) is interpreted by the Lefschetz 
SL(2, C) decomposition of the cohomology of the moduli space M.. They propose the 
Lefschetz decomposition with respect to the fiber SL(2,C)l and the base SL(2, C)r of 
the cohomology H*(M). Once we assume existence of these SL(2, C) actions^ , we may 
arrange this decomposition as 

Ig ® Rg + Ig-l ® Rg-l H + I ® R , (5.5) 

where Ik '■= [(h) ©2(0)]® fc . Identifying this decomposition with that in ( p. If) , they propose 

a k = X (Rk) (k = 0,---,g), (5.6) 

where x(Rk) is the dimensions of the SL(2, C)r representations in R k weighted with 
(— 1) 3R . As argued in [11] we can determine a g and «o easily by the geometric Euler 
numbers; 

a g = (-l) d X (M) , a Q = (-l) d x(M) , (5-7) 

where d = diuuM and d = dinuM. 

Now let us consider our genus g pseudo-sections with a fixed section, say the zero 
section. Counting BPS states for these pseudo-sections are summarized in the generating 
function Q g .\ in ([4.19Q . As we have seen in the last section the genus g pseudo-sections 
come with the moduli M.(g) := Sym s (P 1 ) = P 9 . Furthermore the data of the Jacobian 
may be specified by g points on S, one for each elliptic fiber. Therefore we naturally come 
to the space 

M(g) = S^n 9 (S) , (5.8) 

where ~ represents the resolution of the orbifold singularities via the Hilbert scheme of g 
points on S, which we denote S^K The construction of the Lefschetz actions on this space 
and making the decomposition (|5.5| ) would be interesting problem, however we already 
have predictions for the decomposition; 

I g X x(Rg) + Ig-l X x(Rg-l) + • • ■ + I X X (-Ro) ■ (5-9) 



' One of the Lefschetz SL(2, C) actions is the multiplication of the Kahler form k. Since our 
moduli spaces have a natural fibration structure n : Ai(g) — > M(g) with a section t, we may 
decompose the Kahler class k into fci = (k — ir*(i*(k))) and kn = i*(k). This decomposition 
defines the SL(2, C)i x SL(2, C)r actions. We would like to thank Y. Shimizu for pointing this 
out to us. 
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Namely we can read ( |5.9| ) in our expansion ( |4 . 1 9| ) : 
g = l -2I 1 + 12I 

g = 2 3/ 2 - 30/i + 90/o 

g = 3 -4/ 3 + 52/ 2 - 260/i + 520/ (5.10) 

g = A 5/4-78/3 + 507/ 2 - 1690/i + 2535/o 

g = 5 - 6 h + 108 h - 840 / 3 + 3636 / 2 - 9090 h + 10908 / 

In the next subsection we propose a natural generalization of Gottsche's formula for the 
Poincare polynomials of and reproduce the above predictions. Our generalization of 
Gottsche's formula suffices to determine the decomposition (|5.5|) . 



5.2. Gottsche's formula with SL(2, C) L x SL(2, C) R 

Gottsche's formula describes the generating function for the Poincare polynomials of 
the Hilbert scheme of g points on a surface S. In our case it appears as a natural resolution 

of the symmetric product Sym. 9 (S) for the rational elliptic surface. If we assume the 
existence of the Lefschetz actions SL(2, C)l x SL(2, C) r , then the Poincare polynomial 
can be written 

Pt(SW) = (t 29 ) Tr HHSl!!]) t 2 ^ + ^ (5.11) 
in terms of the diagonal SL(2, C) action. Then the problem is to recover both left and 
right charges in the above formula, namely, Pt L ,t R {S^) = (t 9 L t 9 R )Trt 2 ^ 3 ' L t 2 R 3 ' R . In the case 
of g = 1, we note that the decomposition is unique as follows 

(~ l) LjR © 8(0, 0) L R = (1) L+R © 9(0) L+R (5.12) 
Now let us recall Gottsche's formula[9] 

G(t,q) = T\ ^ . (5.13) 

^ (1 - t 2n ~ 2 q n ) (1 - t 2n q n ) 10 (1 - t 2n + 2 q n ) 

As has been interpreted in [26] there is a close relation between (co) homology elements and 
the bosonic oscillators associated to each elements in H*(S). Under this correspondence 
the classical cohomology is represented by the lowest modes, say ak(— 1) (k = 1, • • • , 12), 
and generate the symmetric product of H*(S). The higher mode excitations afc(— m) come 
from the singular strata of the point configurations. Here it is natural to assume that the 
higher mode excitation ak(—m) have the same spin as the lowest mode a^(— 1), whose spin 
contents are uniquely determined in ( |5.12|) . Under this assumption it is easy to recover 
the SX(2, C)l x SL(2, C)^ spin weights in Gottsche's formula 

G(t L , t R , q) = J] I _ ( tLtR)n - lqn) (1 _ (T^pyj 

" (5-14) 



(1 - t 2 L {t L t R y- l q n ) (1 - tK^)-- V) (1 - (tLt R ) n q n y 
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which provides us the Poincare polynomial Pt L ,t R (S^) as the coefficient of q 9 . Then our 
predictions (5.10) for the Lefschetz decomposition should be verified in the formula 



(t L t 



Rj 



(S [9] ) 



t R =-l 



X(R 9 ) [t L + — +2] +x(R 9 -i) [t L + -+2] +■■ ■+x(Ro) 



(5.15) 

For lower g, we have found complete agreements of our predictions from the generating 
functions Q g -i with those coming from the above formula. For example, for g — 3 we obtain 
from (|5.15| ) 



4 [t L + — + 2 + 52 t L + — + 2 



260 \t L + ^- + 2\ + 520 , 

which should be compared with ( J5.10| ). We may summarize the above results into a formula 
relating Gottsche's formula with our generating functions ^. 18 ): 



G(tL,-l,-£) = ^^;l(? 3 )(tL + ^+2y . 



(5.16) 



5.3. Topological string partition function 

So far we have fixed a section to discuss the moduli space of the curves C g . To recover 
the contributions from the Mordell-Weil lattice, we simply need to multiply the Eg theta 
function to the functions we have discussed. Therefore for the BPS state counting on the 

3 

rational elliptic surface S we consider the function 0£ g (3£, tj)G(— i^, —1, |-) in terms of 
the generalized Gottsche's formula. Now it is easy to deduce the following relation; 

e Ea (3t,t^)G(-t L , -l, f.) =V4i( g )K - i + 2)« 

. (5.17) 

= ^Z 9;l(g )A 2 »- 2 (2 S m|) 2 , 

3>0 

where t^ = e lX and we have used the form of the degenerated instantons given by ( f4.13|) . 

3 

This implies the function 0£ g (3t, tj)G(— —1, provides us the all genus topological 
partition function. We may write ( |5.17| ) explicitly by 



q ^m^~ n (i - t L #»ra - ±w - g z - liq)x {2sm 2 ] ■ (5 - 18) 

Here we recognize the helicity generating function in the left hand side. Thus the topolog- 
ical partition function has a simple but suggestive form that the genus zero function Zq i 
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multiplied by the helicity generating function, which is actually the starting point of the 
analysis done in [11]. 

Here for completeness we prove the compatibility of the above result (|5.17| ) with our 
holomorphic anomaly equation (|4.11| ). For this purpose we note the following identity 



which can be proved in a straightforward way; 



A/2 
sinA/2 



n>l V 



(1 



,n\4 



e iX q n ) 2 (l - e~ iX q n ) 2 



exp 



k>i 



2k 



(5.19) 



The compatibility may be easily verified if we use 3-E2(<? 3 ) = 2^>((/) 2 + E 2 (q) found in 
( fOjip and the value ((2) = since we have |g for the both side of ( |5.18|) after the 
differentiation with respect to E2{q)- Also the formula ( |5.19|) explains the simplification 
we have encountered in Proposition 4.2. Namely we have 



E z *M x29 = Z °M ex P ( 2 E CJ jt e ^ 3 ) (A 

g>0 \ k>l 



2k 



(5.20) 
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Appendix . Picard-Fuchs equations of the mirror X w 

Following [12] [13] the Picard-Fuchs differential operators about the large complex 
structure limit are determined to be 

Pi = 96l - W x 6 y - 66 X 6 Z + 249 z 9 y - 166 2 z - 27x(36 x + y + 26 z + 2)(36 x + 6 y + 26 z + 1) 

- 3y(6 x - 86 z )(6 y - 6 Z ) + z(606 x + 326 z + 32){36 x + 6 y + 26 z + 1) 
V 2 = 6 2 y + y(36 x + 6 y + 26 z + l)(6 y - 6 Z ) 
V 3 = (6 Z - 6 y )6 z - z(39 x + 6 y + 26 z + 2)(36 x + 6 y + 26 z + 1) 

where 6 X = x-^, etc. Looking at the characteristic variety of this system we have deter- 
mined the discriminant ( |3.21| ). 
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